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Accuracy Analysis on Large Blocks of High Resolution Images 
INTRODUCTION 
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“TODAY EXISTING HIGH RESOLUTION 
IMAGES ARE ENTERING INTO 
COMPETITION WITH AERIAL 
PHOTOGRAPHY FOR REGIONAL 
MAPPING PROGRAMS AND OTHER 
EXTENSIVE MAPPING APPLICATIONS 
WHERE HIGH RESOLUTION IS 

REQUIRED” 


OUR STUDY IS CONCENTRATED MAINLY ON 

QUICKBIRD IMAGES 
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Accuracy Analysis on Large Blocks of High Resolution Images 
INTRODUCTION 
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USE OF VERY HIGH RESOLUTION IMAGES 

Mapping, Risk Management, Forest, Geology, Regional Planning, Utlity Corridor 
Planning, Mapping for E911 applications, Defense Mapping, etc. 

1. Information contents 

Rule of thumb: for topographic maps 0.05 - 0.1 mm pixel size in map required 

based on lm pixel size -> map 1:10 000 (~ 1“=800‘) 
based on 0,6m pixel size -> map 1 : 6000 (1“=500‘) 

For orthoimages 8 pixel / mm (0.125 m pixel size) 

based on lm pixel size -> orthoimage 1 : 8000 (~1“=667‘) 
based on 0,6m pixel size -> orthoimage 1 : 4800 (1“=400‘) 

2. Accuracy 

Required accuracy of x, y- coordinates: ~ 0,25mm in map, in orthoimage ~ 0,3mm 

-> required ground accuracy 2,5m / 1,5m (0.25 mm at above map scales) 
Image orientation must be better because of additional error components 
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Accuracy Analysis on Large Blocks of High Resolution Images 
Introduction 
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ACCURACY DEPENDS ON: 

■ Orientation Model used 

" GCPs (Number, Distribution, Geometric & Radiometric Quality, etc.) 

■ For Orthos (DTM accuracy, Resolution, Fidelity of Terrain Reconstruction, etc) 

■ Others 

USED ORIENTATION MODELS IN THE INVESTIGATION: 

• SIMULTANEOUS BUNDLE BLOCK ADJUSTMENT (using ephemeris & 
quaternion algebra) 

• RATIONAL POLINOMIAL COEFFICIENTS (+image correction functions) 

USED IMAGE: 

" A block of 40 QB Basic images with limited overlap (hilly, equatorial, jungle 
zone) 

GROUND CONTROL: 

Manually Transferred from existing oriented aerial photography (1:20,000) 
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Bundle Orientation using ephemeris and quaternion algebra 
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z image 


y image 


x image 
(orbit direction) 



= yo + Xp 


or 


0 = -c 


a u j ( Xp - Xo j ) + a u j ( Yp - Yo j ) + a u J ( Zp - Zo j ) 
a 3l J (Xp - Xo j ) + aj (Yp - Yo j ) + aj (Zp - Zo j ) 


y = 


-c 


a 2l (Xp - Xo j ) + a 22 (Yp - Yo J ) + a 23 J (Zp -Zo j ) 
aj (Xp - Xo j ) + a 32 (Yp - Yo j ) + aj (Zp - Zo j ) 


It is possible to use quaternion multiplication to perform a rotation about an 
arbitrary unit axis m by an angle 0. The quaternion that computes the rotation is: 

q = (a,v ) Where a = cos(©/2) , v=sin(©/2). 

— ^ 

Then q | =|a x sin(0/2), A point in 3D P = (0,P) q=a+bi+cj+dk 


q 2 =Py sin(®/2) 
q 3 = p z sin(©/2) 
q 4 = cos(©/2) 


is rotated by p =n*n*n-i 

* * rotated Q P Q 


The orthogonal rotation matrix corresponding to a 
rotation by the unit quaternion p=a+bi+cj+dk is given 
by q*q- 1 = 


ci + b 


2 c 2 


q*q-i = 2 ad + 2 be 

2 bd - lac 


2 be - lad lac + 2bd 

a 2 -b 2 +c 2 -d 2 2cd — 2ab 
lab + led a 2 -b 2 - c 2 + d 2 


Assuming smooth stable paths (that generally is), then no 


Easting or X e phemeris corrections are applied, only attitude corrections are 
° computed in the adjustment. 


Document number- 6 




Bundle Orientation using ephemeris and quaternion algebra 
The QB Camera Sensor Model includes the Coordinate Systems: 
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Image Coord. System 
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Bundle Orientation using ephemeris and quaternion algebra 
Derivation of the co-linearity equations 
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As the Level IB product is sampled at a constant rate, the corresponding time (t)= t=r/ avgLineRate + firstLineTime 

One point on the imaging ray is the perspective centre of the virtual camera at time t. The coordinates of the perspective centre in the spacecraft 
coordinate systems are constant and given data. In matrix notation: C s = (CX, CY, CZ) T (from *.geo file) 

It is possible to locate the origin of the spacecraft coordinate system in the ECF system at a time t by interpolating the position time series in the 
ephemeris file. Let us call this position S E (t) Likewise, we can find the attitude of the spacecraft coordinate system at a time (t) in the ECF system 
by interpolating the quaternion time series in the attitude file. This quaternion, q s E (t), represents the rotation from the ECF system to the spacecraft 
body system at time t. Then using quaternion algebra, the position of the perspective centre at time t in the ECF coordinate system is: 

C E (t) = (q E s (t))-l C s q E s (t) + S E (t), or 

C E (t) = q s E (t) C s (q s E (t))-l + S E (t) or CE(t) = R E s (t) C s + S E (t) This is the Projection Center if ECEF 

Any point that we measure (c,r) is expressed in the detector system as 

XD=0, YD=-c*detPitch (detPitch: Calibrated distance between pixels) 

In the camera system is X c = cos(a)* XD - sin(a)* YD + Sh x ; 

Y c = sin(a)* XD + cos(a)* YD + Sh Y ; Z c = C (Nominal Principal Distance) 

As Level IB images do not have lens distortion the image point is identical to the measured image point, hence :X C ,=X C , Y C ,=Y C , Z C ,=Z C 

The unit vector wC that is parallel to the external ray in the camera coordinate system is just the position of (XC’, YC’, ZC’) relative to the 
perspective centre at (0, 0, 0), normalized by its length. In matrix notation, this vector is: W c = (X c , , Y c , , Z c , ) T & w c = W c / 1| W c || 

It is possible to convert this vector first to the spacecraft coordinate system and then to the ECF system. The unit quaternion for the attitude of the 
camera coordinate system, i.e., the quaternion for the rotation of spacecraft frame into the camera frame q c s , is in the geometric calibration file 
(*.geo). Then, using quaternion algebra: 

w e = q E s(t) ■' q s c 1 w c q S c q E s(t) or w e = q S e( 0 q C s w c (q S E (0 q C s) '* or usin g matrix algebra, w E = R E s (t) R s c w c 
Hence, the co-linearity equation will be: 

w c = [R E s(t) RS c 1 or w c = (q s E (t) q c s) 1 w e [q s E (t) q c s l T 
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Bundle Orientation using ephemeris and quaternion algebra 
Ephemeris and Quaternion data. Variance-Co-variance Matrix 
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-Ephemeris File: Sample mean and covariance estimates of the position of the 
spacecraft system relative to the ECEF system. These data are produced for a continuous 
image period, e.g., an image strip, and spans the period from at least 4 seconds before 
start of imaging to at least four seconds after the end of imaging 

-Attitude File: Contains Sample mean and covariance estimates of the attitude of the 
spacecraft system relative to the ECEF system. These data are also produced for a 
continuous image period, e.g., an image strip, and spans the period from at least 4 seconds 
before start of imaging to at least four seconds after the end of imaging 

-Sampling rate (timelnterval): Each 0.02 seconds. For intermediate position only linear 
interpolation is required 

-Possibility to implement a weighting schema based on: 

gX 0 2 gX 0 Y 0 gX 0 Z 0 
°Y 0 2 gY 0 Z 0 

°z 0 2 


2 


aqi aq 1 q 2 

aqiq 3 

°qiq 4 

aq 2 2 

aq 2 q 3 

aq 2 q 4 


aq 3 2 

°q 3 q 4 



°q 4 2 
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Block Adjustment using RPC + image correction functions 
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{Image Space (Line, Sample)} FUNCTIONAL RELATION {Object Space (<|>, X, h)} 

p ty- LAT OFF L / - LONG OFF h - HEIGHT_OFF 

LAT SCALE ~ L ON GSCALE HEIGHTSCALE 


y=Line = g((p, X, h) = 


NUM l (P,L,H) 

DEN l (P,L,H) 


a 1 A 

Yb 


NUM s (P,L,H) c T A 
x=Sample- h(<p,X,h)= DEN (p L H) 


NUM t (P, L y H) — u ^ u 2 L u 2 P + a 4 H + u 5 LP + u^LH + u 2 PH + u$L + u^P + u 2 qH + 
a n PLH + 1 i 22 1^ u 22 LP + UjjLH + P + Uj^P^ + u 2 yPH + (i^L H + u^P H + u 2 qH ^ — u A 


DEN , (P, L, H) = b,+b 2 L + b 3 P + b 4 H + b,LP + b t LH + b-PH + b„L 2 + b„P 2 + b,„H 2 + 
b„PLH + b l: L‘ + b n LP 2 + b It LH 2 + b, s L 2 P + h u P' + b,.PH 2 + b u L 2 H + b„P 2 H + b,„H 2 =b T A 


The de-normalized image space coordinates {Line, Sample), can also be computed from: 

Line = y. LINE SCALE + LINE OFF Sample = x. SAMP SCALE + SAMP OFF 

Computed de-normalized image coordinates {Line, Sample ) 

^ Observed image Coordinates (Line, Sample ) 
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Block Adjustment using RPC + image correction functions 
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Linef=^ > (( ?s ,X s ,hJ+d^ ) +V Li Samplef = p 0> (y s ,l s ,h s ) + dp 0> + V Si 

With: Line J and Sampled measured image coordinates of a TP or GCP s with coordinates (cps, Is, hs ) 

£(j)((ps, Is, hs ) and p(j) (<ps, Is, hs ) computed de-normalized image coordinates of the GCP or tie point s 
dc,Q) and dp(j) Correction Functions to account for the difference between measured and computed coordinates 
Nabla Li and Nabla Si are random errors. Or: 

g 0> (q> s ,}. s ,lt s ) =g(y,kh) .Um_SCALE + LmE_OFF P m ( <V S >K>K) =h(<p,X,h) . SAMPLE_SCALE + SAMPLE_OFF 


di; = r 0 + /; . Sampled- 1) . Line dp = s 0 + s s Sample+ s, Line 


(p Li + 8 ^dif + dr? + dr? + ^ d<p s + ^ d? + ^ dh s = v® 

dr w ^ ^ ‘ ^ 


dr? 


dr? 


d<? s 


d _?Li 

dl 


^Li 

dh. 


Vsl + 


d<P ^ds (J> +^-ds (j> +^-dsf + ^-dm + dfs ' 


ds ? 0 dsf ' dsf 

X = (A t PA) 1 (A t PL) 


Sep, 


dX+^-dh=v% 

dl. s dh s * 
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Experimental Tests. Block of 40 QB Basic images 
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Full Control 
Distribution (Case A) 


Perimeter Control and randomly distributed GCPs in the center of 

the Block. (Case B) 
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Experimental Tests. Block of 40 QB Basic images 
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Perimeter control only. 
(Case C) 


Relaxed perimeter control. 
(Case D) 




1 

77 x 

^ 7<1 t 

02 


104 

7 

A 

15 



— 69 

▲ 

— r* 7 

A X 

37 — 

A 


A 

}0 757 


759* ' 92 

A A 


95 




1 

X 

Vif < 

00 

IQ 

1 

ix 



76 

A 

7 

M 

X X 


7 

A 

96 A * 3 * 

A 


7 

A 




7(t (17 

▲ 

07 

05 129 

1 

X 

9 

76 

A 

T 

3 




A 

755 700 

A A 


4 

— f-y 

X X 

A 

X 


153 

A 


A 

09 J28 

1 

X 

17 

7i 

A 

>5 

76 

A 

183 795 

* ▲ 

42 

h 

X 


X 

TO 

A 

TO 

A 

7 

A 

#2 


i 

X 

3 

0 

79 

X 


1 

A 

77 





40 40 

A A 

X 


so - <3 

A A 

3 


140 

A 

i. A 101} 

X 


( 

X 

7 

a 9 

n 

; 

A 

T3 


"W 79 — 

A 

35 ST 

a x 

A 

D 

X 

X 

91 

A 

▲ 

16 

15 


s 

X 

9 

83 

X * 

90 

SO 

1 

X 

15 

” } m 204 
x x A 

. JW 53 

▲ X 

* 

r 

X 


7 

X 

15 7 

A 



5 

1 

X 

11 

7 

5 




. 

X 

■ 

X 

7 

6 

6 

A 

56 

A 




1 

77 A 

* 7 o r 

02 


rod 

1 

A 

15 



— 69 

▲ 

a’-a' 

37 — 

A 


A 

,0 757 


759* * 92 

A A 


95 




7 

A 

7i< 1 

00 

IQ 

7 



76 

A 

7 

— M 

A A 


“7 

A 

96 

A 1 


7 

A 

^757 a 7M 




A 

07 

no fOO 

7 

A 

9 

76 

A 

, 3 

3 

1 

A 

— 76 97- 

A A 

— 3 

A 

A 


753 

A 

X 

7oo /yo 
i A 


A 

09 725 

A 

7 

A 

17 

7i 

A 

>5 

76 

A 

3 753 795 

A A 

AO 41 

" A 

A 


A 

15 

A 

TO 72T 
A 

— 7 
A 

12 


f 

A 

( 

3 

0 

79 

A 


1 

A 

77 



^ ff A'\200 

40 40 

A A 

A 


“50 g 

* A 

3 


140 

A 


6 

A 

7 

89 

72 

A 

73 

W 71 — 

A 

35 ST 

A A 


* 

X 

X 

97 

A 

174 7 

A 

76 


£ 

X 

9 

83 A 

X * 

90 

59 

1 

A 

75 

15 

, m jn 204 

x x A 

. 53 

▲ A 

A 

15 

A 


7 

A 

15 7 

A 


t 

A 

■ 

A 

7 

6 

& 

A 

56 

A 

5 

7 

A 

71 

A 

7 

5 




Control only in the corners of the 
block. (Case E) 
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Experimental Tests. Block of 40 QB Basic images 

Bundle Adj with Ephemeris & Quaternions 











4 


J 

J 

i 



Case A (89) Case B (33) Case C (25) Case D (9) Case E (5) 

■ RMSEp (CGPs) ■ RMSEp (CHK) 


Max Errors on GCPs & CHKPts 










H 

4 


Ll 

4 

1 

j 


Case A (89) Case B (33) Case C (25) Case D (9) Case E (5) 

[WmAX Error on GCPs ■ MAX Error on CHKPts 


Case 

#GCPs 

RMSEp 

#ChKPts 

RMSEp 

A 

89 

0.74 



B 

33 

0.63 

56 

3.13 

C 

25 

0.50 

64 

3.34 

D 

9 

0.47 

80 

5.42 

E 

5 

0.35 

84 

8.71 


Bundle Adjustment with Ephemeris & 
Quaternion. Accuracies in GCPs & Check 

Points 


Case 

#GCPs 

Max Err 

#ChKPts 

Max Err 

A 

89 

3.27 



B 

33 

1.87 

56 

8.61 

C 

25 

1.58 

64 

10.01 

D 

9 

1.45 

80 

19.01 

E 

5 

0.81 

84 

26.49 


Bundle Adjustment with Ephemeris & 
Quaternion. Max Errors on GCPs & Check 

Points 
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Experimental Tests. Block of 40 QB Basic images 

Block Adj with RPCs + Affine image Correction 



Case A (89) Case B (33) Case C (25) Case D (9) Case E (5) 


■ RMSEp @ GCPs 

■ RMSEp @ ChkPts 


Case 

Num 

GCPs 

RMSD [m] 

Max Errors [ml 

Sx 

Sy 

Axmax 

Aymax 

A 

89 

0.21 

0.25 

2.53 

4.93 

B 

33 

0.19 

0.23 

2.51 

4.56 

C 

25 

0.16 

0.22 

2.48 

1.50 

D 

9 

0.16 

0.21 

2.45 

1.48 

E 

5 

0.15 

0.20 

2.06 

1.35 


Internal Accuracy. Statistics on TPs. Bundle Adj 
with Ephemeris & Quaternions 


Case 

#GCPs 

RMSEp 

#ChKPts 

RMSEp 

A 

89 

0.98 



B 

33 

1.07 

56 

2.21 

C 

25 

1.41 

64 

2.93 

D 

9 

1.57 

80 

4.16 

E 

5 

1.70 

84 

5.97 


Block Adjustment with RPCs + Affine Image 

Correction 


Case 

Num 

GCPs 

RMSD [m] 

Max Errors [ml 

Sx 

Sy 

Axmax 

Aymax 

A 

89 

0.29 

0.32 

1.83 

2.57 

B 

33 

0.65 

0.74 

2.01 

3.34 

C 

25 

0.81 

0.88 

2.49 

2.25 

D 

9 

0.88 

0.97 

2.85 

2.98 

E 

5 

0.91 

1.12 

3.08 

3.36 


Internal Accuracy. Statistics on TPs. Block Adj 
with RPCs + Affine Image Correction 
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Analysis of the results 


Accuracy Comparison. Bundle Adjustment using 
ephemeris + quaternion vs. Block Adjustment 
based on RPC and affine image correction 


CASE 

#GCPs 

Ephem. + Quarter. 

RPCs + Affine 

RMSE@ 

GCPs 

RMSE@ 
Chk Pts 

RMSE@ 

GCPs 

RMSE@ 
Chk Pts 

A(89) 

0.74 


0.98 


B(33) 

0.63 

3.13 

1.07 

2.21 

C(25) 

0.50 

3.34 

1.41 

2.93 

D(9) 

0.47 

5.42 

1.57 

4.16 

E(5) 

0.35 

8.71 

1.70 

5.97 


RMS Errors at GCPs and Check Points 


RMSE @ GCPs & Chk Pts [m] 
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Accuaracy Comparison 



— ■ Ephem, + Quarter RMSE @ GCPs 

Ephem. + Quarter RMSE @ Chk Pts 

— ■ RPCs + Affine RMSE @ GCPs 
RPCs + Affine RMSE @ Chk Pts 
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Self-calibration through Additional Parameters 
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an (Xp - Xo j ) + anJ (Yp - Yo j ) + an j (Zp - Zo j ) , , 

= -c — — — — + add par am 

an 1 (Xp - Xo J ) + a 3 2 J (Yp - Yo J ) + a 33 (Zp - Zo 1 ) 


y = 


aj (Xp - Xpi ) + aj (Yp - Yo [ ) + aj (Zp - Zo 7 ) 
aj (Xp - Xo j ) + aj (Yp - Yo j ) + aj (Zp - Zo 7 ) 


+ add par am 



x= 


x= 


y=y+Pl*y y= 

x=x+P2*y x= 

x=x+P3*x*y 
y=y+P4*x*y 
y=y+P5*SIN (y*0.12566) 
y=y+P6*COS (y*0.12566) 
y=y+P7*SIN (y*0.06283) 
y=y+P8*COS (y*0.06283 
y=y+P9*SIN (y*0.03100) 


y+P12*COS (y*0.01600) 
x+P13*SIN (y*0.03100) 
x+P14*COS (y*0.03100) 
x+P15*SIN (y*0.01600) 
x=x+P16*COS (y*0.01600) 
x=x+P17*SIN (x*0.11) 
*SIN (y*0.03) 
x=x+P18*x*y*COS(K) 
y=y+P 1 8*x*y * SIN (K) 


Systematic image errors QuickBird Atlantic City, NJ 
left: overall effect right: without dominating angular affinity 


To avoid over-parameterization the AP are automatically 
eliminated using a stochastical model: 


1. For each AP compute: 



(Tpi — - \j C[ii . (7 o. ti ^ l y 


reject if otherwise 


2. Compute cross-correlation coeffs. for the parameters 


y=y+P10*COS (y*0.03100) y=y+P19*x*y 
y=y+Pll*SIN (y*0.01600) x=x+P20*y*y 


Rij - —j^== ; Rij>0.85 

yjqaqjj 


Then eliminate the parameter 
with smaller value 


ADDITIONAL PARAMETERS IMPLEMENTED IN 3 * Com P ute B = 1 ~ ( dia § N * dia § N' 1 )' 1 ? eliminate the AP 

THE UNIVERSITY OF HANNOVER PROGRAM that B ( > or = 0.85 

BLASPO V 
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Self-calibration through Additional Parameters 
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a. GCPs Manually and automatically transferred from existing 0.5 m GSD Orthophoto 

Height of the GCPs through interpolation on a existing ± 1.3 feet Accuracy DEM (1”=1600’) 
Same DEM for production of Orthophoto 


Digital Ortho. 0.5 m GSD 


QB Basic Image, ~0.6 GSD 
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Self-calibration through Additional Parameters 
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Type of 
Observations 

Number of 
GCPs 

Standard 

Deviation 

[microns] 

Accuracy at GCPs 

Accuracy at Check Pts 

RMSEx 

RMSEy 

RMSEx 

RMSEy 

Manual 

174 

11.2 

0.71 

0.66 



Automatic 

398 

10.1 

0.48 

0.45 



Automatic 

25 

12.3 

0.52 

0.71 

0.69 

0.72 

Automatic 

20 

13.1 

0.59 

0.75 

0.69 

0.88 

Automatic 

15 

13.8 

0.61 

0.78 

0.78 

1.04 


Bundle Orientation with automatic selection and elimination additional parameters. 

Basic QB Image. Area of Atlantic City, NJ 
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CONCLUSIONS 
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-Although high attitude frequencies effects are removed at the time of basic image 
generation, low attitude (Yaw) effects are still present in form of affinity/angular affinity. 
They are effectively removed by additional parameters 

- Bundle block adjustment based on properly weighted ephemeris / attitude quaternions 
(BBABEQ) are not enough to remove the systematic effects. Moreover, due to the narrow 
FOV of the HRSI, position and attitude are highly correlated making almost impossible 
to separate and remove their systematic effects without extending the geometric model 
(Self-Calib.) 

- The systematic effects gets evident on the increase of accuracy (in terms of RMSE at 
GCPs) for loser and relaxed ground control at the expenses of large and strong block 
deformation with large residuals at check points. Systematic errors are more freely 
distributed and their effects propagated all over the block. [No functional model for SE] 

— Block adjustment based on RPCs with systematic image correction functions remove 
significantly the affinity deformation of the basic QB images. Although relaxed ground 
control produces less accurate results (in terms of RMSE on GCPs) than BBABEQ, the 
remaining block deformations are much smaller. Increase of absolute accuracy between 
65 to 80% can be reported. 

— The systematic effects are also noticeable in the internal accuracy of the Block (residuals 
at pass/tie points). These are smaller for relaxed ground control in the BBABEQ and 
opposite in the case of RPCs with affine image correction function. 
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RECOMMENDATIONS 
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Further Studies need to be carried out to be more conclusive: 

1. To conduct accuracy studies on large blocks of images with larger long overlap. 
This will allows a 3D accuracy study and to increase the reliability of the derived 
quantities. 

2. To extend the model of the BBABEQ to include Self-calibration making use of the 
co-variance matrices of the Position and Attitude to build-up proper weighting 
matrix. 

3. To include cross track image strips to study the combined effect of multi-rays 
points and self-calibration on the final accuracy. 

4. To extend the RPC model with other image correction functions, including 
statistical test to ensures the validity of the image correction function parameters 
to avoid over-parameterization. 

5. To test alternative orientation models such as those based on feature matching, 
i.e., matching of linear features existing on a land data base and extracted from 
the image space. 

6. Others 
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